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Abstract 

We study the finite size effects on Bose-Einstein condensation (BEC) of an ideal 
non-relativistic Bose gas in the three-sphere (spatial section of the Einstein 
universe) and in a partially finite box which is infinite in two of the spatial 
directions (infinite slab). Using the framework of grand-canonical statistics, we 
consider the number of particles, the condensate fraction and the specific heat. 
After obtaining asymptotic expansions for large system size, which are valid 
throughout the BEC regime, we describe analytically how the thermodynamic 
limit behaviour is approached. In particular, in the critical region of the BEC 
transition, we express the chemical potential and the specific heat as simple 
explicit functions of the temperature, highlighting the effects of finite size. These 
effects are seen to be different for the two different geometries. We also consider 
the Bose gas in a one-dimensional box, a system which does not possess BEC 
in the sense of a phase transition even in the infinite volume limit. 
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1. Introduction 

The phenomenon of Bose-Einstein condensation (BEC) is described in statis- 
tical mechanics text-books (e.g. [H, [l])- Given an ideal gas of particles obeying 
Bose statistics inside a box with sides of length Li, L2, L3, as the thermody- 
namic limit is taken in the usual way {N -^00, V ^ 00 while N/V and Li/ Lj 
are held fixed), the fraction of particles in any excited state i, N.JN , goes to 
zero. This is expected since the single particle energy levels get closer and closer 
to each other, eventually forming a continuum in the thermodynamic limit. The 
ground state is the exception. Indeed, below a certain critical temperature, the 
fraction of particles in the ground state, Ngr/N, will be non-vanishing in the 
thermodynamic limit, which means that the probability distribution of particles 
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as a function of their energy has a Dirac S component at zero energy. Because 
of this, the chemical potential and, as a consequence thermodynamic functions 
in general, are non- analytical at the critical temperature, signalling the BEC 
phase transition. 

In finite systems, however, we are away from the thermodynamic limit and, 
strictly speaking, a phase transition does not happen, all thermodynamic func- 
tions being smooth functions of the temperature at the critical point. Neverthe- 
less, a large finite system can be practically indistinguishable from an infinite 
one. For example, in the thermodynamic limit as described above, the specific 
heat for the three-dimensional gas has a sharp peak with a discontinuous first 
derivative at the onset of BEC. For a finite large system, there will still be a 
more or less sharp peak but all thermodynamic functions will be analytical. 
The larger the volume that the finite system has, the sharper the peak. The 
difference between a finite system and its infinite counterpart can be brought 
out analytically in the form of finite size correction terms. Furthermore, these 
finite size corrections are dependent on the particular way the gas is confined; 
namely, they depend on the geometry of the system. 

It is our aim in this article to obtain analytically finite size corrections to the 
thermodynamics of an ideal non-relativistic Bose gas in the quantum degenerate 
regime in three distinct physical situations: the three-sphere (spatial section of 
the Einstein universe), the three-dimensional box which is infinite in two of the 
directions and finite in the other one (we will call this the infinite slab) and 
the one-dimensional box. The choice of these models enables us to perform a 
full analytical treatment showing very explicitly in simple expressions the finite 
size effects and the impact that different geometries can have in these effects. 
The three seemingly disparate situations under consideration have mathematical 
aspects in common, allowing for a unified treatment as will be seen. Actually, it 
is well known that the one-dimensional box does not possess BEC, but once we 
have all the mathematical apparatus set up, the study of this system requires 
no extra effort. The infinite slab is an example of generalized BEC, in which the 
particles condense into a low lying set of states rather than the ground state, as 
described early on in Q (and which the work in Q already hinted at) and later 
systematically studied by van den Berg and collaborators dig (see for 
recent articles on this topic). From the systems we study here, the only one 
having the usual form of BEC is the three-sphere. 

Our point of departure in the analytical treatment of the thermodynamic 
sums will be the Mellin-Barnes transform, a tool used in the past in 



I12I [13[ , to study respectively the ideal Bose gas in a harmonic oscillator po- 
tential, in the infinite flat space subject to a magnetic field and in the three 
dimensional space in which one of the directions is compactified to form a cir- 
cle. This transform can also be used to obtain high temperature expansions for 
the ideal Bose gas in quite general settings, as was shown for the Bose gas under 
arbitrary background potentials and in boxes of arbitrary shape 14, 1^ and. 



more recently, for the Bose gas in product manifolds JJ)]. Although these high 
temperature expansions can provide a way of calculating a BEC critical tem- 
perature in each setting, they cannot be used to study the BEC transition itself 
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or the BEC regime. Our approach ahows us to obtain large size expansions, in 
terms of temperature and chemical potential, which are valid and very effective 
throughout the quantum degenerate regime and in the vicinity of the critical 
region. This is done in section For the first two systems mentioned above, we 
will then obtain, in section [31 an analytical description of the approach to the 
critical behaviour which signals the onset of BEC in the thermodynamic limit. 
Specifically, in the critical region we obtain the chemical potential, fraction of 
condensed particles and the specific heat as explicit functions of the tempera- 
ture only. From the results thus obtained, it will be clear that the geometry of 
the system has a crucial role in the nature of the finite size corrections in each 
situation. For example, while in the three-sphere the specific heat peak happens 
at a slightly higher temperature than the temperature at which the chemical 
potential is zero, the opposite happens in the infinite slab scenario. 

The case of the one-dimensional box has been treated before by Pathria 17 1 
who obtained finite size corrections to the number of particles (as a function 
of temperature and chemical potential) using the Poisson summation formula. 
However, his procedure is valid only away from the quantum degenerate regime. 
Finite size corrections to non-relativistic BEC in the three-sphere have also been 
treated before by Altaic [31 • His procedure is similar to that of Pathria's article 
mentioned above. In particular, it is based on the use of the Poisson summation 
formula and likewise, it is valid only away from the quantum degenerate regime. 
A physical system analogous to our infinite slab case was considered before in 
the context of the thermal Casimir effect in 1^, 2^, where the authors obtain 



high temperature (classical region) expansions for the thermodynamic potential 
and tackle the quantum degenerate regime by setting the chemical potential 
identically equal to its (thermodynamic limit) critical value. Some developments 
were also made in the finite size effects of a relativistic Bose gas in the three- 
sphere [U These works used a combination of analytical and numerical 
techniques and were very much inspired by an earlier study by Singh and Pathria 



231. Other work on BEC for the relativistic Bose gas on the 3-sphere includes 
241 . Some work on finite size effects for the uncharged relativistic Bose gas on 



the 3-sphere have been considered 



2. Number of particles and specific heat 

Consider an ideal Bose gas with particle eigenstates of energy En- The 
grand-canonical average number of particles is 

N = ^[e^(^^.-A') _ (1) 

i 

where /3 is the inverse of the temperature T, fi is the chemical potential and the 
sum is over all particle eigenstates. We will use the natural units system with 
h = 1, c — 1, and Boltzmann's constant k = 1 throughout. Expression ^ gives 
us n{T) implicitly if wc fix N. The internal energy is 

C/-^i?,[e^(^--^)-l]-i. (2) 
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The specific heat at constant volume is obtained by differentiating U with re- 
spect to T with the volume and number of particles held fixed: C = [dU /dT)\^ y. 
In this section, we shall aim at obtaining analytical expressions for N and C as 
functions of /i and /3, in the form of expansions that contain the thermodynamic 
limit behaviour plus finite size corrections. 

2.1. The three-sphere 

The energy eigenvalues in the three sphere are given by 

En = - ifi^ + m^a'^ + , n = 1,2,3, ... 

a ^ ' 

with degeneracy g„ = v? (see e.g. [2^ ll^l)- a is the radius of the three- 
sphere, m is the particle mass and ^ is a coupling constant that describes the 
interaction of the field with the scalar curvature. We will limit ourselves to 
conformal coupling, ^ = 1/6 for simplicity. As noted before this is not 
very restrictive since choosing a different coupling is equivalent to changing the 
mass from m to \rr? -f (61^ — l)/a^]^/^. For large a, this change is very small. In 
particular, in the large volume limit, all couplings will give the same results. In 
the non-relativistic limit, E^ can be expanded in powers of v?' jima)'^ yielding 



i?„ = m+— ^ 1-- — +0 — . (3) 

2ma^ \ 4 \ma/ \\ma/ J J 

The constant term, m, can be absorbed into the chemical potential in ([1} and 
([2]), which amounts in fact to a redefinition of the chemical potential (from /i to 
/i — to) in all expressions, with no effects in the physics of the system. Hence, we 
will drop TO from ^ . In addition, we will consider only the first term inside the 
parenthesis, leaving only En — / {2ma'^). Taking this limit requires n <C ma, 
i.e., little occupation of high energy states ("high" when compared with to), 
corresponding to T <C to. For this to be true in the temperature range we 
are most interested in, the critical region, it suffices having p <^ m?, where 
p is the particle density. The relative error involved in this non-relativistic 
approximation is of the order T/toQ 
Wc have then 

n—l n—1 k—1 



^ To see this, we split the sums in |(T]| and l(2)l into a non-relativistic main part running 
from n = 1 to ?i = p, plus a tail (where the non-relativistic approximation is not valid) 
running from n = p -|- 1 to infinity, p must be a number large enough so that the contribution 
of the tail to the whole sum can be neglected and, at the same time, small enough for the 
non-relativistic approximation to be valid in the first part of the sum. By choosing p = 
ma{T/mY = {2x)-^l'^{T jmy-^l'^ , with 1/4 < r < 1/2, it can be shown that for small T/m 
and small x, the relative error in the first part of the sum is of the order T/m and that the 
tail contribution to the error is exponentially suppressed. (See also next footnote.) 
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where we have changed to the dimensionless variables, 



X — ^ ^ „ , e — —2ma^Li . (5) 

(Note that e > — 1 always, for positiveness of the occupation numbers, and 
that e — >■ — 1 implies iV — >■ oo.) x and e can be seen as the dimensionless inverse 
temperature and chemical potential respectively. The number of particles in the 
ground state is A^gi- — [e'^^*^^^ — l] when expressed in terms of dimensionless 
variables. 

In the thermodynamic limit we have a: — >■ 0, which yields the bulk result 

TVbuik = ^U,/,{e-n^-'^' (6) 

when e > and 

when — 1 < e < oH Lia{z) — X^^^i ^"/"" denotes the polylogarithm func- 
tion. In ([7]), the first term on the right hand side is the ground state contri- 
bution (when (e -I- l)x <C 1) and it is of leading order only if (1 -I- e) oc x^^^. 
Since the volume of the three-sphere is given by ^ = 27r^a^, we have Nx'^^'^ = 
(/3/(2TO))^/^27r-^/3, where p is the particle density. Hence, when taking the ther- 
modynamic limit, Nx^^"^ is a constant. Since Li3/2(e~'^^) is a decreasing function 
of ex, defined only when ex > 0, the maximum value that it can take happens 
when e = and it is given by Li3/2(1) = C(3/2). Hence, if the particle density 
exceeds a certain critical value, pc, given by 



2m J ' 4 ^ V2, 

e must go to negative values (it must go close enough to -1) allowing the ground 
state to take the particles that cannot be accommodated in the excited states 
and BEC will happen. Note that this pc is the same as the one for the standard 
case of an infinite cubic box P, 0]- This is natural, since in the thermodynamic 
limit we have a — >■ oo and the three-sphere will resemble more and more an 
infinite flat space. In fact, by writing 2: in ([6]) and ([7]) in terms of volume and 
temperature, we see that these expressions are exactly the same as those of the 
infinite cubic box. 

From the critical density, given as a function of the temperature in ([8]), we 
have the familiar critical temperature given as a function of the particle density 



■^Tho first order relativistic corrections to expressions ^ and (0 for small T/m are re- 
spectively (15v^/32)(T/m)Li5/2(e-^^)a:-3/2 and (150F/32)(T/m)C(5/2)a;-3/2 . (Hence, the 
relative error involved in the non-relativistic approximation is of order T/m, as previously 
stated.) 
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as Tc ~ (27r/r7i)(p/C(3/2))^/'^. Equivalently, the critical value of x is given by 
as can be seen directly from (|6]). 
For e > 0, from ^ we have that in the bulk, ex is a function of the rescaled 
temperature, T /Tc, only (equivalently, Xc/x). More specifically, it is an increas- 
ing function oi T/Tc, as can be clearly seen from (|6]). Hence, away from the 
condensate region (i.e., e > 0), as we approach the thermodynamic limit, we 
have e oc (for fixed T /Tc). In the condensate regime and fixed T /Tc we will 
have e — )■ —1 in such a way that (e + 1) cx ^/x, so that the condensate density 
is not vanishing . 

Of course, these results are only correct in the thermodynamic limit. We 
now obtain finite size corrections to the number of particles in ^ and ([7]) and 
to the specific heat sums, still to be introduced. 

In [Appendix A we define a class of sums, S{i,j), which appear throughout 



the cases we study here, both in the number of particles and in the specific heat. 
For these sums we have obtained asymptotic expansions for small |ea;|. These 
expansions are given in (jA.5|) and (|A.6p . Thus, from (|4]), (jA.ip and (|A.5|) we 
have 

N = Sil, 0) ^ ^x~'/' £ tll^ _ k) (ex)" + X- Vi(l, , 

where /i(a, e) is a function defined in (|A.4|) . The symbol ~ has the usual 
meaning of asymptotic expansion. Specifically, by truncating the infinite sum 
above at a certain order k = n, we know that the error in N (or in any other 
of the sums S{i,j)) will be of order o ((ea;)"). From the earlier discussion based 
on the bulk expressions © and (O, we know that for large N, the condition 
I ex I <C 1 is satisfied throughout the quantum degenerate regime and also for 
temperatures not much higher than Tc- Therefore, for these temperatures, the 
number of particles should be well approximated by only a few terms from the 
sum in k above plus the isolated term. The functions fi{a,e) are studied in 
[Appendix B[ where we find their meromorphic structure and their values at 
particular points. We finally obtain 
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TT-v/e _ 

— — coth(7rveja 



(9) 

where we have truncated the sum at the third term and for the last term 
we require the value of /i(l,e) from (jB.ll) . (|B.4p and (IB.6p . Note that this 
expression is valid for both positive and negative e, although for e < the 
last term is more conveniently written with a regular cotangent by noting that 
•yecoth(7r-y/e) — -y^cot(7r-y^). 

As mentioned in the introduction, the non-relativistic ideal Bose gas in the 
three-sphere has been treated before by Altaie In particular, using a pro- 
cedure which is valid in the e > region only (i.e., away from the quantum 
degenerate regime), Altaie obtained an approximate expression for N which 
consists of the first term inside the square brackets and the term outside the 
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square brackets in © , which amounts to truncating the asymptotic series for N 
above at the fc = term. Hence, our result shows that this is simply the leading 
part of a full asymptotic expansion which is valid for both positive and negative 
e (i.e., also in the BEC regime, where in fact it is most useful) and which can 
be made as accurate as one wishes (at least, in a temperature window around 
Tc) by truncating at an appropriate order. 

When we take the thermodynamic limit we have x — > and we recover the 
continuum approximation given by ^ and ([7]), as would be expected. To see 
this, note that the expansion in square brackets is the small ex expansion of 
Li^/2{e~'^^) obtained by Robinson in except that the term —2^/TT(ex)^^'^, 
present in the polylogarithm expansion, is not present here. However, when 
e — >■ oo (which happens for T > Tc when TV — > oo) the hyperbolic cotangent 
tends to 1 exponentially fast. Hence, the last term in (|9]) will account for the 
missing polylogarithm expansion term and N will be given by ([B]). On the other 
hand, in the BEC regime we will have e ^ — 1 when N oo and, expanding 
y/ecoth^TTy/e) in powers of (e + 1), the last term in © can be written as 



■ coth{^T^/e)x ^ = 



1 3 fn^ 1 , , 

+ -T (£ + 1) + 



+ 1 4 \12 16 

Thus, to first order and when e ~ — 1, is given by ([7]). The -3/4 term in 
this expression gives us the first order finite size correction to the number of 
particles. 

We now turn our attention to the specific heat. In the case of non-relativistic 
particles in the three-sphere, by using ([2]), ([5]) and the energy levels En = 
/ (2ma^) with degeneracy gn = n^, we have for the internal energy 



1 °° —1 



n=l 

For the specific heat at constant radius a, we obtain 

s{2,iy 



C^x 



5(1,1) 



(11) 



where we used [dN/dx)\^ ^ = to obtain {d{xe) / dx)\^ ^. Inside the square 
brackets we have the omnipresent sums S{i, j) defined in (jA.ip . which once again 
have asymptotic expansions for small ex given in (jA.sp . Altaic also obtained 
expressions for the specific heat sums, but these are only valid in the e > 
region as they yield an imaginary specific heat when e < 0. 

In figure 1, we show the fraction of particles in the ground state, A^gr/TV, as 
a function of the rescaled temperature, T /T^. In figure 2 we show the specific 
heat per particle. For both plots we used © (with only the terms explicitly 
displayed) to obtain the values of e to input in A^g,. (for figure 1) and in the 
sums that appear in (jlip (for figure 2). For these sums, we used their asymptotic 
expansions (IA.5I) truncating the k sums ai k — 2. In both figures, we can clearly 
see the approach to the familiar thermodynamic limit critical behaviour as TV 
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is increased. In the next section, we provide an analytical description of this 
phenomenon. Note that contrary to what happens in the case of a Bose gas 
in a harmonic trap, in our case the finite size effects bring the features which 
are characteristic of the BEC transition to higher values of the temperature. 



This also happens in the case of a Bose gas in a box [151, |29| or in a power-law 
trap with power higher than 3 [s^l. Had we also plotted the exact numerical 
results given by (|4]), ([TT]) and (lA.ll) the figures would not look any different 
because the numerical plots would totally superimpose the analytical ones. In 
fact, our analytical approximation is extremely good for temperatures near 
and throughout most of the condensate region. Taking as an example the N = 
10^ case, the relative error in the condensate fraction, (7Vgi./Af(analytical) — 
iVgi./A^(numerical))/(A^gr/A''(numerical)), is never more than 3 x 10~^ for 0.15 < 
T/Tc < 1.3 being of the order of 10~^ dXT — Tc- For temperatures lower than 
Q.lbTc it deteriorates due to the increase in a; (x — oo as T — )• 0), while 
e ~ — 1. As the temperature rises above Tc, the approximation deteriorates due 
to the increase in ex. At T = 2Tc or T = O.OlTc the relative error is already 
6 X 10"^ (still quite smaU). For TV = 10^ the error is similar in the classical 
regime (since ex is largely insensitive to N in this regime) but much smaller in 
the condensate regime (of the order of 10""'^° or less in the critical region and 
in most of the condensate regime). The relative error in the specific heat has a 
similar behaviour, although it is somewhat larger (very roughly, by a factor of 
about 10 in most regions). Even higher precision is easily attained by truncating 
the asymptotic series at a higher order. Less precision is attained if we truncate 
at a lower order. Had we truncated the asymptotic expansion for iV at fc = 0, 
we would get a relative error in the condensate fraction of order 0.1 already for 
T = 1.2Tc, independently of N . In the condensate region, it would be a few 
percent (order 10~^) in the N = 10^ case and of order 10~^ in the N — 10^ 
case. 

2.2. The one- dimensional box 

We now consider an ideal Bose gas in a one-dimensional box with Dirichlet 
boundary conditions. In this case, the non-relativistic energy levels are 



TT 



where L is the length of the box. The levels are non-degenerate. After the 
change of variables, 

e = -^M, 



and using ([ij, we have 



N = 



00 _^ 
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Looking at this expression, the only thing that distinguishes this gas from the 
one in the three-sphere is the levels degeneracy. 

In the thermodynamic limit, we have the bulk expression 



bulk 



= ^:z:-V2Lii/2(e— ), 



(13) 



which is valid to first order. In this case, we do not need to take special care 
about the ground state because, unlike the polylogarithm of index 3/2 of the 
previous case, Lii/2(e~'^^) is unbounded when ex — ?> 0. This means that, inde- 
pendently of the value we set the temperature and density at (or equivalently, 
the value we set Nx^^"^ at), all particles can be accommodated in a smooth 
distribution over the energy levels, without any condensation. BEC does not 
happen. Another thing that can be seen from ()13|) is that, in the bulk, ex is 
a well defined function of Nx^^"^ (or in more physical variables, a function of 
p/\/T), similarly to the three-sphere case. In the present case, however, this 
holds for all temperatures. 

Using (|A.6|I for S'(0,0) in ([T2|) we have the asymptotic expansion for small 



N 



'1/2 



24 



-ex + X 



(14) 



where we have truncated both infinite summations in (jA.6|) ai k — 2. 

A very similar expression, which is equivalent to this one, was obtained before 
by Pathria [l7| . making use of the Poisson summation formula in a procedure 
that is valid only for e > 0. Once again, our method extends a previously found 
formula to the e < region. 

Similarly to the first case studied, to first order, this expression just gives 
the bulk result p3)) . Indeed, the expansion inside the first square brackets is the 
small ex expansion of Lii/2(e~'^'^) except for the missing term •y/7r(ea:)~-^/^ (which 
is the term that provides the correct divergent behaviour of the polylogarithm 
when ex — > 0). This missing term is accounted for by the first term inside the 
second square brackets (coth(7r-ye) tends to 1 exponentially fast when e — )• oo). 
Hence, we recover ([13]). The other terms in (IT4l) are of second order, i.e., they 
are finite size corrections. 

Because (|14p is asymptotic for small ex and ex increases with temperature, 
our analytical results for the one-dimensional box are valid only in the low 
temperature (or high density) regime. The range of validity is increased by 
including terms of higher order in ex. 

For the specific heat, we can easily derive 



C 



5(2,1) 



5(0,1) 



(15) 
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The behaviour of the specific heat and ground state population is uninteresting 
in the ID box case since all quantities are always smooth. For this reason, we 
do not include any plots here. 



2.3. The infinite slab 

In this section, we consider a Bose gas contained in the space limited by two 
infinite parallel planes. By this, we mean a gas contained in a rectangular box 
with sides of length Li, L2 and L3, in the limit of L2,L^ — ^ 00, while Li and 
L2/L3 are held fixed. 

The energy eigenvalues of a non-relativistic particle in such a box, with 
Dirichlet boundary conditions, are 



En — 



2m 



Li U2, 



Inserting (|T6l) in ([T]), we have 



N = 



00 00 



2. 

A;— 1 ni 



-kp- 



, ni = 1,2,3. 



(16) 



(17) 



Now, in the limit ^2,^3 — ^ 00 the sums in n2 and ns become integrals^ The 
number of particles will also be infinite in this limit. It is then convenient to 
define a new quantity, 77, which will be useful later, as 77 = N L\/ {L2L^). 7] is 
the number of particles in a cube with sides of length Li. It is finite (as long as 
Li is finite) and we prefer using it instead of the particle density, as 77 is more 
suitable as a mathematical analogue of N in the previous cases. 

Performing the integrals in n2 and 713 and changing to the dimensionless 
variables 

/37r2 



2mLl 



277lL? ' 



we obtain 



(18) 



Again, one of the sums S{i,j) defined in (jA.ip makes its appearance. 

The thermodynamic limit is achieved by taking Li — >■ 00 (with density and 
temperature kept constant). In this limit, we merely replace the sum in n by 
an integral, obtaining the bulk quantity 



r3/2 



^bulk 



Li3/2(e-^^)x-3/^ 



(19) 



^Strictly speaking, the limit operations L2, L-j — > 00 must be outside all sums in 1171 1. That 
they can be performed inside the k and ni sums is assured by the uniform convergence of the 
summand in this limit. 
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Similarly to the previous cases, for finite Li this expression gives us the leading 
behaviour for small x (large Li) when e > 0. The situation is identical to the one 
in the three-sphere, as can be seen by comparing (IT^ with The maximum 
value that the right hand side can take is ( (3/2) (7r/(4x))3/2, when e = 0. As 
we take the thermodynamic limit, rya;'^/^ = p(/37r2/(2m))^/2 is a constant. If this 
constant is set at a value higher than {-k /^^Y^'^C, (3/2), then e will have to go to 
negative values and close enough to -1. In this case, the n = 1 part of the sum in 
P^ . which we denote by 771 and is given by 771 = — (7r/4)x^^ In [l — e"'^^^"'"'^^] , 
has to be taken into account separately, just like in the three-sphere case. This 
yields in the thermodynamic limit, 

« = "^^-"'C (0 - ln(. + 1) , (20) 

when e is close to -1. The last term is of first order if ln(e + 1) cx — a;"^/^. Note 
that, unlike in S^, this term is not the ground state contribution. Instead, it 
represents the 771 contribution to 77 (i.e., the contribution of the particles that 
are not excited in the Li direction). 

If we set the temperature at a value lower than the critical value, (which 
is the same as the one in the three-sphere or the one in the more standard cubic 
box case), the fraction of particles with 7^1 = 1, 771/77, will be non- vanishing in the 
thermodynamic limit. What happens is a generalized BEC The usual 

form of BEC does not happen at any temperature or density. Indeed, the ground 

state occupation number is given by A^g,. = e^^^"^^i/^2+i'i/^3+'^) — 1 . From 

here, we see that if Li is finite and e > — 1 (which is the case in all regimes if Li 
is finite), we will have A^gi- finite, even when L2 and L3 are infinite. Since, when 
L2 and L3 are infinite, N is also infinite, we have Ng^/N = in all temperature 
regimes. This is valid when Li is finite and also, naturally, when Li — 00. This 
last scenario can be viewed as an extreme particular case of the anisotropic 
boxes studied in the thermodynamic limit by van den Berg in 

Inserting the results of [Appendix A| for S'(0,— 1) in (fT8|) . we obtain the 
asymptotic expansion for r; 



■■I, ■■ 1 , 2 sinh(7r-i/e) 

e+-x-'\nx x^'\n ■ (21) 

16 8 4 ^ ' 

Like in the previous two cases, by taking the thermodynamic limit in this ex- 
pression we recover the first order results ((T9)) and ((20)) . 

In what concerns the specific heat, the analogue of expressions (fTT|) and (|15p 
here is 



, . , ^ 2 , , [i^(0.-l) + 5(l,0)]' 
x5(2, 0) + 25(1,-1) + -5(0, -2) - 



(22) 
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Hence, using the asymptotic expansions of [Appendix A yet again yields an 



asymptotic expansion for the specific heat. However, in order to achieve this, 
in the specific case of S'(l, —1) and S{0, —2) we use (|B.9p (see last paragraph of 
[Appendix A[ ), which is valid only for |e| < 1. 

We can now obtain e for given rj from (PT|) and use it to calculate 771 , which, 
as we have seen, is the analogue of the three-sphere Ng,-. The results for rji/r] 
are plotted in figure 1 in the cases 77 = 10^ and r] = 10^. Like in S^, if fully 
numerical results were plotted, they would totally superimpose the analytical 
ones (in fact, the error is even lower in the present case). In figure 2 we plot 
C/N obtained from together with the expansions for the sums S{i,j). The 
point where the curves change the pattern in figure 2 is where e becomes larger 
than 1, rendering our results for S'(l, — 1) and ^(O, —2) not valid any more and 
making us resort to fully numerical calculations (when e > 1). In both figures 1 
and 2, the curves are similar to those pertaining to the three-sphere. In the large 
limit they tend to be the same, specifically, the standard text book curves of 
the infinite box. However, we see that the finite size effects have some differences 
from one situation to the other. Namely, the condensate fraction approaches 
the bulk limit curve faster in S'^ with increasing TV than in the infinite slab with 
increasing rj and the specific heat peak is higher and at a higher temperature in 
the infinite slab than in S^. In the next section, these differences are brought 
out analytically. 



3. Critical region temperature expansions 

Expressions (|9]) and ([2T|) give us e implicitly as a function oi N or rj and from 
here we have all thermodynamic quantities. However, it would be much more 
convenient to have e given explicitly. That is what we aim for in this section, 
concentrating on the critical region. We solve ([9]) and ([21]) perturbatively for e. 
This enables us to obtain expansions for e (and hence, the chemical potential) 
and specific heat in powers of (T/Tq — 1) (where Tq is defined below). The effects 
of finite size and geometry are contained in the coefficients of these expansions. 
In this way, we obtain an analytical description of how the system approaches 
critical behaviour as the thermodynamic limit is approached. 

3.1. The three-sphere 

Define xo and respective temperature To as being the values of x and T at 
which the chemical potential assumes the critical value fj, = (hence, at which 
e = 0). From we have 

Naturally, in the thermodynamic limit we have xq /xc — > 1 and for large particle 
numbers xq will be very close to Xc. More precisely, Tc/Tq = xq/xc ~ 1 — 
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0.292:, 



1/2 



Expanding e in powers of [x/xoY/'^ 

1/2 



1 around a; = we have 



ai 



1 



a2 



1/2 



+ ■ 



(24) 



for some coefficients ai and 02- Using ([9]) we find for ai 

,-1/2 (3V^/4)C(3/2)-xy^ 



ai 



^26 + (^/4)^(l/2)4/2 
For large A^, a;o ^ 1 ^-nd we have to the leading order, 



9C(3/2) „i/2 
- 2^3/2 ^0 • 



(25) 



The result for 02 is quite long and there is no point in displaying it in full. For 
Xo 1 it becomes, to leading order. 



27C(3/2)^ _i 

a-) ~ x„ . 

^ 207r ° 

Inserting (pst and (pS)) in (IM)) . we finally have the approximation 



-2.11X, 



-1/2 







1/2 

— I -1 

Xq 



2.93x, 







1/2 

— I -1 



H 2 



(26) 



(27) 



To make this expansion more readily interpreted, it can be put in terms of num- 
ber of particles and temperature if we use 
in powers of {T/Tq — 1), yielding 



and expand ^Jx/xq = ^JT\JT 



.^0 



Jo 



(28) 



Hence, in a neighborhood of Tq, we have e (or the chemical potential, if we use 
([5|l) given explicitly as a function of the temperature, for a given number of 
particles. In this expression, we clearly see e becoming increasingly steeper at 
Tq as the thermodynamic limit is approached (N — >■ 00), a fact that could be an- 
ticipated from the bulk expressions. Because of the way in which we expanded 
the hyperbolic cotangent in ([9]), these results are valid only for |e| < 1 (and 
of course, for small T/Tq — 1). As can be seen from p8|) . this is satisfied if 
]\['^/^(T/To ~ 1) is sufficiently small. Therefore, the temperature window for 
which (|28p is valid is increasingly narrow as we approach the thermodynamic 
limit. 

We can now find an expansion for the specific heat by inserting p4| in (jA.5[) . 
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which in turn is inserted in (jlip . We obtain 



TV 



4c(i) ^ 

729C(i: 



81C 



45C 



AOtt 

-1/2 

1.93 + 0.69 1 ^-1 

Jo 



14OO0F 



4C(|) 

1/2 



1.25iVi/3 ( 1 

To 



T 



1/2 



(29) 



The coefficient of each term is given only to highest order in Xq (or N^^'^). 
The first term is the famihar specific heat maximum of three-dimensional infinite 
space. The other terms are specific to the three-sphere. From the coefficients 
in ([29l) we can calculate the first and second derivatives of the specific heat 
with respect to T at T = Tq. The first derivative does not depend on N, 
but the second one goes to infinity as — > oo. This is not surprising; it is 
the renowned specific heat peak growing sharper. By differentiating (j29l) with 
respect to T and equating the result to zero, we obtain the position of the 
specific heat maximum to be Tmax/7o ~ 1 + 0.28N^^^'^ , to first order in the 
number of particles. Replacing T^ax in ([28|) . we obtain roughly e ~ 0.3. So, 
^max falls within the region where the approximation is valid. This is confirmed 
by comparison of our analytical results with those of our numerical calculations. 
We also present C/N at Tq expanded for small xq (large N): 



15C(|) 
4C(|) 



15C 



27C(i) 



1/2 



1.93 - 0.797V"^/^ 



(30) 



a formula that also agrees very well with the numerical results. The sub-leading 
term gives us the rate at which the specific heat at Tq is decreasing as we get 
away from the thermodynamic limit. 

By inserting (P7)) or (^5]) in the expression for the ground state population, 
we immediately have the condensate fraction as a function of temperature in 
the critical region. For purposes of comparison with the infinite slab geometry, 
we limit ourselves to presenting the leading behaviour at T = Tq' 



_gr 

N 



V^C(3/2) 



-2/3 



N-^/^ ~ 0.917V"^/^ 



3.2. The infinite slab 

For the infinite slab a similar analysis to that described in the previous 
section can be followed. In this case, xq is still defined as the value of x at 
which fi — and is given from (|2ip as 



r3/2 



-C 



-3/2 



Cq ^ In Xq 



ln(27r) 
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We now have 



0.14x1/2 In 



Naturally, as before lim xq/xc = 1. Comparing this expression for T^/Tq with 

N—^oo 

the same one for the three-sphere, we see that in both cases To is higher than 
Tc- However, due to the presence of the log factor in the present case, we have 
that Tq/Tc is higher in this case than in S^. This is consistent with the fact 
that the specific heat peaks in the infinite slab are to the right of the respective 
peaks in the three-sphere. 
The result for e is 



9C(f) 
7r3/2 



.-1/2 



T 



1/2 



27C 



l-73r;i/3 ( - - 1 



0.9877' 



2/3 



IOtt 

T 



The coefficients are only displayed to highest order in Xq 
heat we obtain 



1/2 



1/2 ] ^ 

-1 +••• 

(31) 

For the specific 



N 



15C 



4C(i) 
243Ci 



81C(|)' 45C(f) 



20n 



-1/2 



4C 



1/2 



— -1 



1400^7^ 



1.93-1.51 



1/2 

-I -1 



T 



0.36/7 



1/3 



T 



(32) 



Comparing p2p with (1291) we see that the qualitative behaviour is very similar 
except that now the second term has a negative coefficient, meaning that Tq 
is to the right of the specific heat maximum, in contrast to what happens in 
the three-sphere. The second derivative of the specific heat with respect to T 
still becomes infinitely negative as ry — )■ oo. Once more, this corresponds to the 
sharpening of the specific heat peak. 

If we try to get the value, Tmax, where the specific heat has its maximum, 
as we did in S'^, we get Tlnax/To ~ 1 — 2.117]^^/"^. However, this result is 
not trustworthy because the value i^^^^^iT^iax/To — 1) ~ —2.11 is quite large, 
rendering ([31]) not valid for T — Tmax- 

The expansion of the specific heat at Tq for small xq is 



N 



15C(i) 15C(|) 



4C(i) 40FC(i)' 
1.93 + 0.34?7"i/^ In 7/ 



■^y^ In 



xo 



(2^)2 
3.0777^1/^ + 



^3/2 



27C 



3C 



47r3/2 



1/2 



(33) 



which is the analogue of ((30|) for the infinite slab. This expansion converges 
more slowly than its equivalent in S'^. For this reason, the expression in p3p is 
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not accurate for low particle numbers like rj = 10^. Naturally, it gains precision 
as 77 increases or as more terms are included. We see that here the leading 
finite size correction is positive for large enough particle numbers (specifically, 
for Inxo ^ —5.6 which corresponds to 77 > 8.6 x 10'^). This is opposite to what 
happens in the previous case, where the leading finite size correction is always 
negative, and consistent with the fact that the specific heat peaks are higher in 
the infinite slab than in S^. 

For the condensate fraction, at the temperature Tq we have to first order 
r]i/rj ~ (2/3)C (3/2)~^^^ ?7~^/^ In?] ~ 0.3577^^/'^ In 77. Comparing with the cor- 
responding expression in the three-sphere, we observe that there is an extra 
logarithmic factor. This means that in the present case the condensate fraction 
at To does not tend to zero as quickly as in S'^ when the thermodynamic limit 
is approached. Indeed, this can be observed in figure 1. 

4. Conclusion 

We studied the finite size effects on the thermodynamics of an ideal Bose gas 
at low temperature, especially in the BEC regime. For this, we obtained asymp- 
totic expansions for the number of particles and specific heat sums which contain 
the bulk behaviour in their leading contributions, while the other contributions 
contain the finite size corrections. The three cases studied have mathemati- 
cal similarities that enabled a unified treatment. The asymptotic expansions 
proved very accurate in the BEC regime and also in the vicinity of the critical 
point, even when keeping only two or three terms. They are not as useful in the 
one-dimensional box due to the absence of BEC in this case. 

We were able to provide expressions for the chemical potential and specific 
heat as explicit functions of the temperature only in the vicinity of the critical 
point, in the form of expansions whose coefficients depend only on the size and 
geometry of the system. These expressions provide an analytical description 
of how thermodynamic functions approach their bulk forms as the systems ap- 
proach the thermodynamic limit, showing for example the gradual sharpening 
of the specific heat peak. Furthermore, they allow us to see clearly the effects 
of finite size and the differences between these effects in different geometries. 
For instance, while in the three-sphere the specific heat maximum happens at 
a temperature T > To the opposite happens in the infinite slab and the leading 
finite size correction to the specific heat value at Tq is negative in the former 
case and positive in the latter. 

It should be possible to apply the methods employed here to study BEC 
in other systems. In particular, it is likely that the relativistic ideal Bose gas 
can be treated in this way. As mentioned in the introduction, finite size effects 
for relativistic bosons in the three-sphere have been studied by several authors 
using a combination of analytical and numerical techniques, but fully 
analytical results have not yet been obtained for that system to the accuracy 
that we have described here. It would also be interesting to consider the finite 
size systems studied here with an interacting Bose gas and see to what extent 
our ideal gas results would be affected. 
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Appendix A. The sums S{i,j) 



In this appendix, we consider the class of sums 

oo oo 

S'(i,j) = ^^n2*/fcJe-'=^("'+'^) (A.l) 

n=l fe=l 

which appear in the expressions for the number of particles and specific heat 
in the different cases studied. We will obtain asymptotic expansions for small 
ex. For our applications, we need to consider the cases i — 0, 1,2, 3 and j = 
-2,-1,0,1. 

Apply the Mellin-Barnes transform 

6-"= = —/ dar(a)a;-", x > 0, c> (A.2) 

to the exponential in (jA.ip . Then, we can move the summations inside the 
integral sign as long as there is uniform convergence, which is the case if we 
choose c > max{0, j + 1, i + 1/2}. The sum in k becomes a Riemann C- function 
(which we denote by C,) and we have 

S{i,j) = 7r- daf,{a,e)T{a)C{a-j)x-", c> max i 0, j + 1, i + 



(A.3) 

where we defined the functions fi by 

oo ^ 

/,(a,e) = ^n2'(„2+e)-", 3?(a) > » + _ . (A.4) 

n=l 

The functions /j are dealt with in [Appendix B[ There, it is shown that they 
can be analytically continued to the whole complex plane except for isolated 
singularities at a = 1/2 + i — k, k — 0,1,2,.... The other factors are also 
analytical in C except for poles at a = j + 1 (coming from ^(a — j)) and 
at a = 0,-1,-2,... (coming from T{a)). We now close the vertical path of 
integration in the integral in (|A.3p in a large rectangle to the left and solve it 
using the residue theorem. For this, we must know the residues of the integrand 
In what concerns fi{a,e), all the information required is given in [Appendix B 



The residues of T{a) at a = —k are (— l)*/fc! and the residue of C,{a) at a = 1 
is 1. Due to the decaying properties of the integrand for large values of 9(a), 
the expansions thereby obtained are asymptotic for small \ex\. These are, for 
i = 1, 2, 3 and j = 0, 1, 

Si^, j) - /,(j + l, + V2 g tllr Q + ( + ^ - J - fc) 

(A.5) 
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For i = and j = 0, 1, we have 

k=0 ■ ^ 

+ 2 5I^^C(-J-A:)(ex)'=. (A.6) 

fc=0 

The first term in (jA.5l) conies from the ^-function pole, while the infinite sum- 
mation comes from the poles of fi. In (jA.6[) it is exactly the same except that 
we have a second infinite summation. This comes from the F-function poles. It 
is not present in (|A.5p because these poles are removable in this case, due to 
the zeros of fi{a, e) when i ^ 0. 

When j = —1,-2 the situation is only slightly different in that the pole 
coming from the ^-function is double and so the term with k ~ — j — 1 is omitted 
from the second summation in (IA.6[) and the first term in (jA.Sp and (jA.6|) will be 
different. Specifically, it will be given by the residue of fi{a, e)F(a)C(a — j)x~°' 
at a = j + 1. This is obtained in a straightforward fashion from the results of 
Appendix B[ more precisely, from (|B.ip . (|B.6I) . (|B.7p and (|B.9p . 



Appendix B. The functions fi{a.,e) 

In [Appendix A[ equation (|A.4p . we have defined the functions fi{a, e). Here 
we obtain their meromorphic structure in the variable a as well as simple ex- 
pressions for fi{a, e) and {d/da)fi{a, e) at some particular values of a. We need 
this information for the calculations in [Appendix A[ 

First notice that fi{a,e) for i = 1,2,3 is easily obtained from /o(q:, e) by 
using the recurrence relation 

f,+i{a,e) = f^{a - l,e) - e/j(a,e). (B.l) 

Hence, we only need to study /o(a, e). Equation (|A.4p defines fo{a,e) as an 
analytical function of a for 3?(a) > 1/2. There is more than one way of obtaining 



the analytical continuation of this function. Ghika and Visinescu [31| , Ford [32] 
and Elizalde and Romeo all found procedures to do it, but only for e > 
(whereas we have e > —1). (Ford's result is actually for a more general function 
of which ours is a particular case.) The same function as in [s^l was treated by 
Actor [s^j, who obtained the analytical continuation in the case |e| < 1. The 
procedure followed in [isj for a similar function would yield the required pole 
structure (in the whole range of e). Here we show a different route. 

The analytical continuation can be achieved by noting that /o(a, 0) = C(2q!) 
and that 

^/o(a,e) = -a/o(a + l,e), (B.2) 

or in its integral version, 

/o(a, e) - C(2a) ~a f /o(a + 1, e) de . (B.3) 
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This immediately provides the analytical continuation to 3fJ(a) > —1/2. In this 
region, we see that the only pole will be the one of C(2q!) at a = 1/2, its residue 
being 1/2. The other poles and residues follow easily by using (|B.3p repeatedly, 
replacing /o(a, e) and /o(a + l,e) by their Laurent expansions. Specifically, 
/o(q;, e) is analytical in C except for simple poles at a = 1/2 — fc, fc = 0, 1, 2, . . . 
and the respective residues are given by 



Resfc(e) ^ Q ^ J Resfc_i(e)de, 



fc = 1,2,.... 



where ReSfc(e) denotes the residue of /o(a, e) at a = 1/2 — fc. By induction, we 
finally obtain 



Resfc(e) = (-ir- --1 •••U-fc 



2 V2 / V2 J kl 2 r(l/2-/fc)fc! 

From (jB.ip it is easily seen that fi{a, e) is analytical in C except for simple poles 
at a = 1/2 + j - fc. A; = 0, 1, 2, . . . 

Expressions for fo{a, e) at a = 2, 1, 0, — 1, . . . are also needed. The case a = 1 
is given by [s^ 

oo 

/o(l, e) = ^(n^ + e)-' = ^ coth(vrVi) - ^ • (B.4) 

n— 1 ^ 

/o(0,e) follows directly from (|R3| with a = yielding /o(0,e) = -1/2. The 
other values are obtained from these ones by using (jB.2p or (|B.3p as needed, 
yielding 

/o(2, e) = ^csch2(^yi) + ^ coth(7rVi) - ^ , (B.5) 

/o(-/c,e) = -ie^ fc = 0,l,2,... (B.6) 

Now, we want the derivatives dfo/da at a = and a — —1. Inserting the 
Taylor series of fo{ct, e) and /o(l + a, e) around a = in (jB.3p . yields 



^/o(a,e) 



^_^^2sinh(Vf)^ (B.7) 



a=0 V 

where we used (|B.4p . Following the same procedure for dfo/da at a = —1, we 
obtain 



^fo{a,e] 



. 1 f\ 2sinh(7rVi) , 
= 2C'(-2) + -e- / In de , (B. 



where the prime in C'(— 2) means derivative. There is no simple form for the 
integral in (jB.Sp . However, in the case of |e| < 1, it can be expressed in terms 
of an infinite sum (see, e.g., [s^), which is more convenient than just leaving it 
as it is. We then have 



^_C(3)_ 2sinh(7r7i) ^ (2^)^fei32fc ^+1 
2^2 ^^^(2 + 2fc)(2fc)!' 



(valid for |e| < 1), where B2k are the Bernoulli numbers. 
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Figure B.l: Condensate fraction Ng^/N in for = 10^ (continuous line) and AT = 10® 

(alternating short/long-dashcd line) and r]i/r} in the infinite slab for r) = 10^ (long-dashed 
line) and r] = 10^ (short-dashed line) as functions of the rescaled temperature. 




Figure B.2: The specific heat per particle in 
(short-dashed line) and in the infinite slab 
long-dashed -numerical- line) and rj = 10^ 
-numerical- line) - see text for details. 



S3 for N = 10^ (long-dashcd line) and JV = lO"' 
for rj = 10^ (continuous -analytical- and then 
(continuous -analytical- and then short-dashed 
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